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Abstract 

Follicle-Stimulating Hormone Receptor {FSHR) -29G/A polymorphism (rsl 394205) was reported to modulate gene expression 
and reproductive parameters In women, but data In men Is limited. We aimed to bring evidence to the effect of FSHR -29G/ 
A variants In men. In Baltic young male cohort (n = 982; Estonians, Latvians, Lithuanians; aged 20.2±2.0 years), the FSHR -29 
A-allele was significantly associated with higher serum FSH (linear regression: effect 0.27 lU/L; P = 0.0019, resistant to 
Bonferroni correction for multiple testing) and showed a non-slgnlficant trend for association with higher LH (0.19 lU/L) and 
total testosterone (0.93 nmol/L), but reduced Inhlbin B (-7.84 pg/mL) and total testes volume (effect -1.00 mL). Next, we 
extended the study and tested the effect of FSHR gene haplotypes determined by the allelic combination of FSHR -29G/A 
and a well-studied variant c.2039 A/G (Asn680Ser, exon 10). Among the FSHR -29A/2039G haplotype carriers (A-Ser; 
haplotype-based linear regression), this genetic effect was enhanced for FSH (effect 0.40 lU/L), Inhlbin B (—16.57 pg/mL) 
and total testes volume (—2.34 mL). Finally, we estimated the total contribution of three known FSH-actlon modulating 
SNPs {FSHB -21 IGAT; FSHR -29G/A, c.2039 A/G) to phenotypic variance in reproductive parameters among young men. The 
major FSH-actlon modulating SNPs explained together 2.3%, 1.4%, 1.0 and 1.1% of the measured variance in serum FSH, 
Inhlbin B, testosterone and total testes volume, respectively. In contrast to the young male cohort, neither FSHR -29G/A nor 
FSHR haplotypes appeared to systematically modulate the reproductive physiology of oligozoospermic Idiopathic Infertile 
patients (n = 641, Estonians; aged 3 1.5 ±6.0 years). In summary, this Is the first study showing the significant effect of FSHR - 
29G/A on male serum FSH level. To account for the genetic effect of known common polymorphisms modulating FSH- 
actlon, we suggest haplotype-based analysis of FSHR SNPs {FSHR -29G/A, c.2039 A/G) In combination with FSHB -211G/T 
testing. 
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introduction 

Follicle-stimulating hormone (FSH) secreted by anterior pitui- 
tary together with other endocrine factors plays a central role in 
establishing and maintaining human fertility. Circulating FSH 
stimulates gametogenesis and steroidogenesis in gonads by binding 
into its receptor (FSHR). During male fetal, neonatal and pubertal 
periods, FSH stimulates proliferation of testicular Sertoli cells 
determining spermatogenic capacity of adult testes, and in 
adulthood it contributes to normal spermatogenesis and sper- 
matogonia! survival and sperm release [1,2]. Inactivating muta- 
tions in the FSH P-subunit coding FSHB and the FSH receptor 
coding FSHR genes result in severely impaired spermatogenesis 
[3,4]. In addition to loss-of-function variants, common polymor- 
phisms in these genes have been shown to contribute to male 



reproductive physiology [5]. We have previously shown that the 
T-allele of the FSHB -21 IG/T promoter variant (rsl 0835638) was 
associated with significantly reduced serum FSH levels and total 
testes volume in the Baltic cohort of young men [6,7], and these 
results were confirmed in Estonian, Italian and German infertile 
male patients [8-10]. Recent studies have also conclusively shown 
the association between the Ser680-allele of FSHR 2039A>G 
(p.Asn680Ser, rs6166) and higher serum FSH, lower total testes 
volume, Inhibin B and total testosterone levels [11,12]. The FSHR 
c.2039A>G and the linked variant FSHR c.919A>G 
(p.Thr307Ala, rs6 165) jointly determine the two FSHR isoforms 
[13]. 

Another common polymorphism FSHR -29G/A (rsl 394205) in 
the 5 '-untranslated region of the gene (Figure lA) has been 
reported to affect its transcriptional activity [14]. Although the 
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discovery study in women liad concluded tliat tliis SNP lias little 
impact on FSHR expression and gonadal function [15], the 
follow-up clinical studies in female patients have showed that the 
expression level of the FSHR on human granulosa cells obtained 
from AA-genotype carriers was only 40% compared to the GG- 
genotype carriers [16] and fittingly, the amount of exogenous FSH 
required for ovulation induction was 1.8-fold higher in AA- 
homozygotes [16,17]. The data on the effect of FSHR -29G/A in 
men is limited. Despite studies on several populations [18-22] and 
meta-analyses across individual reports [20,23] have failed to 
identify significant contribution of the RSHR -29G/A on male 
infertility, there is missing data on the effect of this SNP on a wider 
range of male reproductive parameters. A pilot study has reported 
smaller testicular volume in Estonian men carrying the FSHR -29 
A-allele [20]. 

We set forward to bring conclusive evidence to the effect of 
FSHR -29G/A in men in analysing a large study group (n = 1,623) 
comprising of Baltic young male cohort (n = 982) in comparison 
with Estonian oligozoospermic idiopathic infertile male patients 
(n = 641). There is emerging data on the importance of SNP-SNP 
and gene-gene interactions within relevant biological pathway(s) in 
determination of the studied phenotypic variation [24]. We 
performed haplotype-based association analyses combining the 
FSHR -29G/A genotype data reported in this study with the 
previously published dataset of the FSHR c.2039A>G (p.Asn680- 
Ser) genotypes on the same samples [12]. In addition, the study 
groups were subjected to analysis of the joint contribution to 
normal phenotypic variance of the three main genetic variants 
reported to modulate the FSH action {FSHB -211G/T, FSHR - 
29G/A, FSHR c.2039A>G) [5]. To our knowledge, this is the first 
study showing the significant effect ot FSHR -29G/A alone and in 
combination with FSHR c.2039A>G and FSHB -211G/T on 
male serum FSH level and downstream reproductive parameters. 

Materials and Methods 

Ethics statement 

The study has been approved by the Ethics Committee of 
Human Research of the University Clinic of Tartu, Estonia 
(approval date 27.01.2003), the Ethics Committee of Riga Stradins 
University, Latvia (23.04.2003), and the Regional Ethics Com- 
mittee of Kaunas, Lithuania (approval no. 13, 2003). 

The Baltic young male cohort 

The Baltic male cohort was recruited between May 2003 and 
June 2004 among the participants in a prospective study 
Environment and Reproductive Health (EU 5th FP project 
QLRT-200 1-029 11) in parallel at tiiree study centres (Tartu, 
Estonia; Riga, Latvia; Kaunas, Lithuania). The recruitment and 
phenotyping protocols at the participating centres were identical. 
Study participation was voluntary and written informed consent 
was obtained from all subjects. Details of the study group 
formation were described previously [25]. Men were recruited to 
the study at the Centre of Andrology, University Clinic of Tartu, 
Estonia (n = 578; all born and living in Estonia), at the Riga Family 
and Sexual Problems Centre, Latvia (n = 300; ail born and living 
in Latvia), and at the specialized laboratory' of the Institute of 
Endocrinology, Kaunas University of Medicine (n = 326; all born 
and living in Lithuania). Previously, either the sub-cohort of 
Estonian men or the full cohort of Baltic young men has been 
investigated for two SNPs modulating FSH-action, FSHB -21 IG/ 
T [6,7] and FSHR c.2039A>G [12]. In genetic association studies 
we have excluded cohort participants with clinical factors leading 
strongly deviated reproductive physiology (lack of sperm in 



ejaculate, i.e. azoospermia, n = 2; cryptorchidism, n=13; abuse 
of anabolic steroids, n= 1; orchitis with unilateral testis damage, 
n=l) or incomplete clinical data (n = 15). In addition, for the 
current study the DNA samples of a subset (n = 190) of the fuU 
Baltic male cohort were not available. The final number of 
analysed Baltic young male cohort participants for FSHR -29G/A 
(rs 1394205) was 982 (Table 1). 

Estonian idiopathic infertility patients 

The study group of oligozoospermic Estonian men with 
idiopathic infertility (n = 750) was recruited at the Andrology 
Centre, Tartu University Clinics between June 2003 and August 
2008 and consisted of male partners of couples failing to conceive 
a child for a period of S12 months. Oligozoospermia was 
diagnosed according to the World Health Organization (WHO) 
criteria valid at the time of recruitment (sperm concentration 
<20 mln/mL [26]). Phenotyping protocol was identical with that 
in Baltic young male cohort [25]; the details of the formation of 
the study group are described elsewhere [27]. In brief, all study 
participants were of white European ancestry, born and living in 
Estonia. All men with causal factors for male factor infertility 
(obstruction, cryptorchidism, chromosomal abnormalities, Y 
chromosome deletions, hypogonadotrophic hypogonadism, testic- 
ular diseases, sexual dysfunctions, androgen abuse, severe traumas 
and operation in genital area, chemo- and radiotherapy) were 
excluded from the analyses resulting in a study group consisting of 
688 participants. For the current project, also patients with 
azoospermia, i.e. lack of sperms in ejaculate, n = 47) were 
additionally excluded from the genetic analysis. The final number 
of analysed subjects successfully genotyped for FSHR -29G/A 
(rs 1394205) was 641. Among the included infertility patients 
(sperm concentration <20 mln/mL), 408 subjects also fulfilled the 
latest WHO 2010 criteria for oligozoospermia (sperm count below 
39xl0Vejaculate [28]). 

Hormone assays 

For all participants of the study, venous blood was obtained 
from the cubital vein in the morning and serum was separated 
immediately. Serum sampling period for the Baltic cohort was 
from 08.00 to 13.00 (median 1 1.00), and for the Estonian infertility 
patients from 08.00 to 11.00 (median 9.30), respectively. 

For the Baltic cohort, serum levels of FSH, LH and total 
testosterone were determined using time-resolved immunofluoro- 
metric assays (Delfia, Wallac, Turku, Finland), estradiol by 
radioimmunoassay (Pantex, Santa Monica, CA, USA) and Inhibin 
B by a specific two-sided enzyme immunometric assay (Serotec, 
Oxford, UK) at the Department of Growth and Reproduction in 
Copenhagen, Rigshospitalet, Denmark in the framework of the 
Environment and Reproductive Health (EU 5th FP project 
QLRT-2001-02911). The intra- and inter-assay coefficients of 
variation (CV) for measurement of both FSH and LH were 3 and 
4.5%, for total testosterone <8% and <5%, for estradiol 7.5% 
and 13%, and for Inhibin B 15% and 18%, respectively. 

For the Estonian idiopathic infertility patients, the FSH, LH, 
total testosterone and estradiol levels of blood serum were 
measured using the Immulite automated chemUuminescence 
immunoassay analyser (Immulite; Diagnostic Products Corp., 
Los Angeles, CA) according to manufacturer's instructions, at the 
United Laboratories, University of Tartu Clinics. At the same 
laboratory, Inhibin B was determined in duplicate using a specific 
enzyme immunometric assay (Diagnostic Systems Laboratories, 
Inc., Webster, TX). The intra- and inter-assay CV were 4.2 and 
8% for FSH; 4.0 and 7.1% for LH; 6.3 and 9.4% for testosterone; 
7.5% and 13% for estradiol; 15% and 18% for Inhibin B. 
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Figure 1. Genomic structure of the FSHR and effect of the FSH-action modulating genetic variants on reproductive parameters. A, 

Schematic representation of the structure of the FSHR gene drawn to an approximate scale. Exons are depicted as boxes, translated sequences are 
shaded in grey and transcription start-site is denoted with up-right arrow. Circle-headed bars indicate the location of FSHR SNPs -29G/A and 4-2039A/ 
G (Asn680Ser); and their mutual genomic distance and linkage disequilibrium is provided. B, Distribution of the four FSHR haplotypes (G-Asn, G-Ser, 
A-Asn, A-Ser) formed from the FSHR -29G/A and +2039A/G (p.Asn680Ser) variants in the Baltic male cohort and Estonian oligozoospermic infertility 
patients (number of phased chromosomes, n = 1964 and n = 1282, respectively). C, Individual allelic effects of the FSHR -29 A-allele and FSHR 680Ser 
(c. 4-2039 G), and the cumulative haplotypic effect of the formed FSHR A-Ser gene variant on serum FSH, Inhibin B and total testes volume in the Baltic 
young men cohort. Results of the association testing are presented as P-values and effect sizes (regression coefficient, P) from linear regression 
analysis. Arrows indicate the strength and direction of the effects. D, Proportion of total phenotypic variance (%) of serum FSH, Inhibin B, total 
testosterone and total testes volume explained by the FSHB -21 IGAF (white bars), FSHR -29G/A (grey bars) and FSHR Asn680Ser (black bars) genetic 
variants in the Baltic young male cohort. Individual and cumulative effects of the SNPs were estimated by using the REML analysis implemented in 
GCTA software [33], 
doi:1 0.1 371 /journal.pone.0094244.g001 
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Table 1. General characteristics of the study groups. 





Parameter 


Baltic male cohort (n = 982) 


Estonian oligozoospermic patients (n = e41) 


General characteristics 


mean±SD 


median(5-95th percentile) 


mean±SD 


median(5-95th percentile) 


Age (years) 


20.2 ±2.0 


19.8(17.4-24.2) 


31.5±6.0 


30.9(23.4-42.0) 


BMI (kg/m^) 


22.3 ±2.5 


22.1(18.8-27.0) 


26.6±4.4'' 


25.9(20.7-34.9) 


Abstinence period (hours) 


107.7±63.2 


86.0(48.0-231.4) 


92.3±51.0 


72.0(48.0-168.0) 


Total testes volume (mL) 


49.2±10.3 


50.0(33.0-70.0) 


40.3 ±10.3 


40.0(24.0-56.0) 


Sperm concentration (10^/mL) 


81.7±74.4 


63.3(9.2-214.3) 


7.8±5.9 


7.0(0.1-18.0) 


FSHR -29G/A (rsl 394205)'" 


Allele frequencies 


G 


74.8 (1466) 




77.1 (988) 




A 


25.4 (498) 




22.9 (294) 




X^-test' 






P = 0.12 




Genotype frequencies 


G/G 


56.2 (552)" 




59.3 (380)' 




G/A 


36.9 (362) 




35.6 (228) 




A/A 


6.9 (68) 




5.1 (33) 




X^-test= 






P = 0.12 





'Data for BMI available for 324 patients of the oligozoospermic study group. 
''Data presented as percentage with number of allele/genotype carriers in brackets. 

'^P-value from x^'l^^^t for differences in FSHR -29G/A allele and genotype distribution between Estonian oligozoospermic infertile patients and Baltic male cohort. 
"Hardy-Weinberg Equilibrium test P-value of =0.4. 
'Hardy-Weinberg Equilibrium test P-value of =1.0. 
doi:1 0.1 371 /journal.pone.0094244.t001 



Testing of genetic association between tlie FSHR -29G/A 
(rs 1394295) SNP and serum hormone levels was carried out 
separately for the Baltic young men cohort and Estonian 
idiopathic infertility patients. 

Semen analysis and physical examination 

Semen samples were obtained by masturbation and all semen 
values were determined in accordance with the World Health 
Organization (WHO) criteria valid at the time of recruitment [26]. 
In brief, after ejaculation, the semen was incubated at 37°C for 
30-40 min for liquefaction. Semen volume was estimated by 
weighing the collection tube with the semen sample and 
subsequently subtracting the predetermined weight of the empty 
tube assuming 1 g = 1 mL. For assessment of the sperm concen- 
tration, the samples were diluted in a solution of 0.6 mol/L 
NaHC03 and 0.4% (v/v) formaldehyde in distilled water. The 
sperm concentration was assessed using the improved Neubauer 
haemocytometers . 

Patients were examined by clinical investigators who had passed 
special clinical training. Physical examination for the assessment of 
genital pathology and testicular size was performed with the man 
in standing position. If necessary, pathologies were clarified further 
with the men in supine position. The orchidometer (made of birch 
wood, Pharmacia & Upjohn, Denmark) was used for the 
assessment of testicular size. The total testes volume is the sum 
of right and left testicles. 

Genotyping procedure and data 

Genomic DNA was extracted from peripheral blood using a 
modified version of the salting-out method [29]. The FSHR -29G/ 
A (rs 1394205) polymorphism was genotyped by PGR and allelic 
discrimination assay on the ABI PRISM 7900HT detection system 
(Applied Biosystems, Foster Gity, GA, USA). Each PGR reaction 



(5 |il) contained 20 ng genomic DNA, 1 |il 5x HOT FIREPol 
Probe qPGR Mix Plus (ROX) (Soils BioDyne, Tartu, Estonia), 
0.25 |il predesigned TaqMan SNP Genotyping Assay (ID: 
G_426553_10; Applied Biosystems, Foster Gity, GA, USA) and 
3.75 |il ddH20. Default thermal cycling conditions were applied: 
15 min at 95°G followed by 40 cycles of 15 sec at 95°G, 1 min at 
60°G, and 2 min at 50°G. Results were analysed using the allelic 
discrimination endpoint analysis mode of the Sequence Detection 
software package, Version 2.4 (SDS 2.4). Quality control of 
genotyping was guaranteed by inclusion of positive control DNAs 
(previously sequenced) representing the three FSHR -29G/A 
genotypes (GG, GA, AA) into each genotyping plate. 

Previously published genotyping data for FSHR Asn680Ser 
(c.2139 A/G, rs6166; Table SI; [12]) and the FSHB -211G/T 
(rs 10835638; Table S2; [7]) were integrated into the FSHR 
haplotype (-29G/A-C.2139 A/G) analysis and in estimating the 
cumulative effect of all the three FSH-action modulating SNPs on 
phenotypic variance of the studied parameters. 

Data analysis 

Mean, standard deviation, median, and 5-95* percentiles were 
calculated for general characteristics (age, BMI, ejaculation 
abstinence period) and main outcome variables (hormonal and 
seminal variables, total testes volume) using PASW software Grad 
Pack 18.0 (SPSS Inc., Ghicago, IE, USA). Genepop software 
(Version 4.0.10) was used test concordance of the genotyping data 
with Hardy-Weinberg equilibrium, and to test population 
diflFerentiation (x2-test) between the two study groups [30]. The 
testing conditions were: dememorization = 10000, batches = 1000, 
iterations = 10000. 

Genetic associations between the studied individual SNPs and 
male reproductive parameters were tested using multiple linear 
regression implemented in PEINK, version 1.07 (http://pngu. 
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mgh.harvard.edu/purcell/plmk/) [31]. The natural log-transfor- 
mation was used to obtain an approximate normal distribution of 
values for all studied quantitative parameters except total testes 
volume. Single SNP a.ssociation tests were performed under 
additive genetic model. In the Baltic male cohort, regression testing 
was performed with the adjustment for age, BMI, smoking status, 
and recruitment centre. Hormone measurements were additionally 
corrected for blood sampling hour, and semen parameters were 
corrected for ejaculation abstinence period according to the analysis 
settings described previously [7]. In case of Estonian idiopathic 
infertile patients, linear regression was performed with adjustment 
for age; semen parameters were additionally corrected for 
abstinence period. The Bonferroni threshold for correction for 
multiple testing in linear regression analysis was estimated 0.05/ 
16 = 3.13x10"^, taking into account the number of independent 
measurements (eight) and tested study samples (two). 

The amount of LD between the FSHR -29G/A (rs 1394205) and 
Asn680Ser (rs6166) was calculated using r" and D' statistics 
implemented in Haploview [32]. The FSHR gene haplotypes 
comprised of alternative alleles of the FSHR -29G/A and c.2039 
A/G (p.Asn680Ser) and subsequent hapl()t\pc frcqu(mcics were 
inferred from unphased genotype data using exjxx tation-maximi- 
zation (EM) algorithm implemented in PLINK, ver. 1.07 [31]. 
Tests for associations of the inferred FSHR -29G/A - c.2039 A/G 
haplotypes and male reproductive traits were conducted using 
haplotype-based linear regression analysis implemented in PLINK 
by ' — hap-hnear' option based on a sliding-window approach with 
two included SNPs. Haplotype-based analyses were adjusted for 
appropriate cofactors analogously to the individual SNP associa- 
tion tests. Correction for multiple testing was performed using 
Max (T) permutation procedure (permutations = 10,000). Haplo- 
type omnibus tests were performed testing for an overall effect of 
the haplotypes on the studied parameter. 

Non-parametric Mann- Whitney two-tailed [/-test implemented 
in PASW software Grad Pack 18.0 (SPSS Inc., Chicago, IL, USA) 
was used in Baltic male cohort to assess statistical difiFerences in 
serum FSH level and total testes volume between the carriers 
(GA-f-AA, n = 430) and non-carriers (GG homozygotes, n = 552) 
and of the FSHR -29 A-allele. Mann-Whitney J7-test compares the 
medians and the distribution of values. The Bonferroni threshold 
for correction for multiple testing in Mann-Whitney f/-test was 
estimated 0.05/8 = 6.25 X 10 '\ taking into account the number of 
independent measurements (eight). 

The proportion of total phenotvpic variance of studied 
hormonal and testicular parameters explained by the FSHB - 
211G/T, FSHR -29G/A and FSHR AsnGSOSer SNPs was 
estimated using the REML (restricted maximum likelihood) 
analysis implemented in Genome-wide Complex Trait Analysis 
(GCTA) software (http://www.complextraitgenomics.com/ 
software/gcta/index.html) [33]. Briefly, Genetic relationship 
matrixes (GRMs) were calculated for each study sample to 
determine the genetic relationship between pairs of individuals. 
GRMs served as input into a restricted maximum likelihood 
(REML) analysis to produce estimates of the proportion of 
phenotypic variance explained by the studied SNPs (Vg/Vp). 
REML analyses were adjusted for appropriate cofactors analo- 
gously to the individual SNP association tests. 

Results 

FSHR -29 A-allele is not enriched among male infertility 

patients 

FSHR -29G/A SNP (rs 1394205) was genotyped in the Baltic 
male cohort (n = 982; age 20.2 ±2.0 years; sperm concentration 



8 1.7 ±74.4 mln/mL) and in Estonian oligozoospermic patients 
diagnosed with idiopathic infertility (n = 641; age 31.5±6.0 years; 
sperm concentration 7.8±5.9 mln/mL; Table 1). The estimated 
allele and genotype frequencies of the FSHR -29G/A did not differ 
between the two study groups (minor allele frequency, 
MAP = 25.4% in Baltic cohort vrs. 22.9% in Estonian oHgozoo- 
spermic patients; Fisher's exact test, P=0. 12; Table 1). 

FSHR -29G/A is a novel identified genetic determinant of 
serum FSH 

Among the Baltic young men, linear regression analysis resulted 
in significant association of the A-allele of the FSHR -29G/A SNP 
with higher serum FSH (additive model; P= 0.0019, A-allele effect 
0.27 lU/L, resistant to Bonferroni correction; Table 2). This was 
supported by statistical analysis comparing the two subgroups of 
patients stratified based on their carrier status of the FSHB -29 
A-allele. The group of consisting of AA-homozygotes and GA- 
heterozygotes (n = 430) compared to GG-homozygotes (n = 552) 
had significantiy higher serum FSH (median 2.9 vrs. 2.6 lU/L, 
Mann- Whitney {/-test P= 0.004, resistant to Bonferroni correc- 
tion; Table 3). Additionally, the A-allele of the FSHR -29G/A 
showed a non-significant trend (P<0.06) for association with 
increased serum LH (linear regression: effect 0.19 lU/L) and total 
testosterone (effect 0.93 nmol/L) levels, but lower Inhibin B 
(effect —7.84 pg/mL) and total testes volume (effect —1.00 mL) 
(Table 2). Serum estradiol and semen parameters (semen volume, 
sperm concentration and count) were not associated with the 
FSHR -29G/A variant. 

In contrast to healthy young male cohort, no genetic 
associations between FSHR -29G/A and tested reproductive 
parameters reached Bonferroni-corrected statistical significance 
level irrespective whether oligozoospermia was diagnosed accord- 
ing WHO 1999 criteria used at the recruitment (sperm concen- 
tration below 20xlO''/mL [26], Table 2) or based on the revised 
WHO 2010 edition (sperm count below 39xl0'^/ejaculate [28], 
Table S3). Interestingly, among infertilit}' patients, the FSHR -29 
A-allele showed a trend for association with higher total testes 
volume (oligozoospermia based on WHO 1999, nominal 
P= 0.017, Table 2; WHO 2010, P= 0.043, Table S3). 

Gene haplotypes formed from FSHR -29G/A and c.2039 
A/G (p.Asn680Ser) exhibit enhanced effect on serum 
hormones and testes volume compared to individual 
SNPs 

We aimed at the haplotype-based association analyses combin- 
ing the FSHR -29G/A genotype data reported in the current study 
with our previously published tiie FSHR c.2039A>G (p.Asn680- 
Ser) genotype data of the same samples [12] (Table SI). The 
common polymorphisms FSHR -29G/A (in 5'UTR) and FSHR 
c.2039 A/G '(p.Asn680Ser, exon 10) kicate >190kb apart and 
showed linkage disequilibrium (LD; D' = 0.027; = 0.0) neither in 
the Baltic male cohort nor in Estonian oligozoospermic patient 
study group (Figure lA). The four FSHR common variants arise 
from alternative allelic combinations of the FSHR -29G/A and the 
FSHR c.2039 A/G SNPs: -29A/2039A (A-Asn), -29A/2039G 
(A-Ser), -29G/2039A (G-Asn), and -29G/2039G (G-Ser). The 
estimated frerjuency distribution of the four haplotypes did not 
differ between the young men and oligozoospermic male patients 
(Table 4, Figure IB). The most prevalent haplotype was -29G/ 
2039A (G-Asn) with 44.6% and 45.9% carrier frequency in the 
two study groups, respectively. 

We addressed the overaU effect of the FSHR -29G/A - c.2039 
A/G inferred haplotypes (allelic combinations of these two 
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Table 2. Marker-trait association analysis and clinical parameters of the two study samples stratified into subgroups based on the 



FSHR -29G/A (rsl 394205) genotypes of the participants. 




Parameter 




Baltic male cohort* 






Estonian oligozoospermic men* 








Group 


mean±SD 


median(5-95%) 


/'-value 


beta(SE)'' 


mean±SD 


medlan(5-95%) 


/'-value 


beta(SE)'' 


FSH (lU/L) 


G/G 


3.0±1.6 


2.6(1.1-5.9) 


0.0019* 


0.27(0.07) 


7.7±6.6 


5.6(1.9-21.6) 


0.098 


-0.41(0.26) 




G/A 


3.2±1.8 


2.9(1.2-6.6) 






6.9±5.3 


5.5(1.9-17.5) 








A/A 


3.4±1.7 


3.0(1.4-7.0) 






6.3 ±4.5 


4.7(1.2-17.8) 






1 H Ml I/I ) 


G/G 


3.9±1.7 


3.7(1.8-6.9) 


0 021 


n 1 9(0 081 


4.5±2.2 


4.1(1.6-8.5) 


0 50 


— 0 ORfn 1 7) 




G/A 


4.2 ±1.6 


4.0(1.9-7.1) 






4.2 ±1.8 


3.9(1.8-7.5) 








A/A 


4.1 ±1.5 


3.7(2.1-7.0) 






4.5 ±1.9 


4.0(1.7-8.3) 






Inhibin B (pQ/mL) 


G/G 


234.6±83.1 


225.0(115.0-391.6) 


0.039 


— 7.84 (3.90) 


87.3±54.9 


78.2(17.5-179.4) 


0.093 


10.60(6.12) 




G/A 


224.6 ±7 1.9 


219.0(110.2-351.9) 






100.9±58.6 


92.0(11.5-197.0) 








A/A 


217.3±72.1 


198.0(109.7-363.2) 






11 7.5 ±106.7 


78.0(15.6-374.5) 






Total testosterone 
(nmol/L) 


G/G 


27.1 ±9.5 


25.9(14.5-44.3) 


0.042 


0.93(0.45) 


18.6±5.9 


18.0(10.5-29.0) 


0.40 


-0.35(0.42) 




G/A 


27.8±8.4 


26.9(14.6-42.9) 






18.6±7.3 


17.9(8.8-31.9) 








A/A 


28.6 ±9.9 


27.5(16.3-51.1) 






18.2±5.3 


17.7(10.1-30.5) 






FttraHini lr\mr\\/\ \ 

a\J\\-}\ \IJI \ \\Jl/ l-J 


G/G 


92.9±25.8 


89.0(58.0-141.0) 


0 26 


1 35(1 20) 


99.2±36.1 


87.0(73.0-155.8) 


0 55 


1 21 (2 01) 




G/A 


95.0±23.9 


92.0(59.3-142.7) 






101.8±40.2 


88.1(73.0-171.3) 








A/A 


94.6±26.3 


90.0(60.0-135.3) 






100.0±31.9 


90.2(73.0-1 70.0) 






Total testes volume 


G/G 


49.5 ±10.9 


50.0(32.0-70.0) 


0.057 


-1.00(0.53) 


39.4±9.6 


40.0(22.0-55.0) 


0.017 


1.67(0.69) 




G/A 


49.2 ±9.3 


50.0(36.0-65.0) 






41.5±11.4 


41.0(25.0-59.3) 








A/A 


46.2±10.4 


47.0(27.2-62.2) 






41.6±10.1 


41 .0(24.2-62.4) 






Semen volume (ml) 


G/G 


3.5±1.6 


3.4(1.3-6.4) 


0.55 


-0.05(0.08) 


4.2±1.7 


3.9(1.7-7.7) 


0.53 


0.07(0.12) 




G/A 


3.6±1.7 


3.3(1.3-6.6) 






4.2 ±1.8 


4.0(1.7-7.7) 








A/A 


3.4±1.7 


3.0(1.4-6.5) 






4.8 ±2.3 


4.9(1.6-9.7) 






Sperm concentration 
(loVmL) 


G/G 


85.5±82.2 


64.3(9.2-225.4) 


0.28 


-2.94(2.85) 


7.8±6.0 


7.0(0.1-18.0) 


0.37 


0.43(0.47) 




G/A 


77.4±65.4 


63.0(8.2-202.0) 






7.5±5.7 


7.0(0.1-18.0) 








A/A 


72.5 ±48.0 


62.6(7.1-148.9) 






10.4±6.5 


10.4(0.1-19.0) 






Total sperm count 
(10'^/ejaculate) 


G/G 


288.5 ±298.3 


219.5(19.1-780.8) 


0.21 


-11.79(9.91) 


33.1±31.2 


24.5(0.4-95.0) 


0.30 


1.99(1.91) 




G/A 


264.0±252.4 


198.3(19.1-719.0) 






32.7±29.6 


23.0(0.3-89.4) 








A/A 


235.3± 184.5 


195.9(18.1-662.0) 






48.3±37.6 


45.6(0.3-126.3) 







"Baltic young men cohort, n = 982, A-allele frequency 25.4%, HWE test P = 0.40; Estonian oligozoospermic men, n = 641, A-allele frequency 22.9%, HWE test P=1.0. 
^FSHR -29 A-allele effect is shown as the estimated linear regression (additive model) statistic beta (P), standard error of the regression (SE) is shown in brackets. Asterisk 
(*) indicates a significant association, P<0.05 after Bonferroni correction for multiple testing. 
doi:1 0.1 371 /journal.pone.0094244.t002 



polymorphisms) on male hormonal and testicular parameters, as 
well as analysed the individual effect of the four haplotypes 
(Table 4). Overall, in both study groups the carrier status of FSHR 
haplotypes exhibited a significant effect on serum FSH (haplotype 
omnibus test: Baltic cohort, P= 0.007; Estonian oligozoospermic 
infertile patients, P= 0.045), and total testes volume (P<0.03). In 
young men, the FSHR haplotypes also affected significantly serum 
Inhibin B levels {P= 0.043). Among individual haplotypes, the 
least common FSHR variant (A-Ser; frequency 9.8%) in young 
men was associated with significantly higher serum FSH 
(P= 0.0033, haplotype effect size 0.40 lU/L), but lower Inhibin 
B (P= 0.022, effect -16.57 pg/mL) and total testes volume 
(P= 0.021, effect -2.34 mL). The joint effect of the two genetic 
variants on each parameter was cumulative, e.g. FSHR -29 A allele 
effect on FSH is 0.27 lU/L, FSHR Ser-variant effect 0.06 lU/L 
and their haplotype effect 0.4 lU/L (Figure IC). In case of total 



testes volume, the effect of minor alleles of these SNPs and their 
formed haplotype is additively negative: — 1 .00, —1.13 and 
— 2.34 mL, respectively (Figure IC). In contrast, the most 
prevalent FSHR G-Asn gene variant (44.6%) was positively 
correlated with larger total testes volume (P=:0.015, effect 
1.34 mL). 

The results of inferred haplotype analysis are in line with the 
data of young men stratified based on the carrier status of nine 
alternative genotype combinations of the FSHR -29G/A and 
Asn680Ser polymorphisms (Figure 2). The carriers of both, FSHR 
-29A and 680Ser-variants exhibited significantiy higher level of 
serum FSH (Mann- Whitney Z7-test P= 2.5 x 10"'^; Figure 2A) and 
a trend for lower total testes volume (P<0.1; Figure 2B) 
compared to the rest of the study subjects. The homozygotes for 
the minor alleles of both FSHR variants (AA-SerSer) were 
measured 13% lower total testes volume compared to the subjects 
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Table 3. Test results for the difference in estimated median values and distributions of the study parameters in the Baltic male 
cohort subgroups stratified based on the carrier status of the A-allele of FSHR -29G/A. 



Parameter* 


FSHR -29 GG-carriers (n = 552) 


FSHR -29 A-carrlers (GA+AA, n = 430) 


/'-value" 


FSH (lU/L) 


2.6(1.1-5.9) 


2.9(1.3-6.6) 


0.004* 


LH (lU/L) 


3.7(1.8-6.9) 


3.9(2.0-7.1) 


0.016 


Inhibin B (pg/mL) 


225.0(115.0-391.6) 


218.0(110.6-352.5) 


0.091 


Total testosterone (nmol/L) 


25.9(14.5-44.3) 


27.0(14.9-45.0) 


0.049 


Estradiol (pmol/L) 


89.0(58.0-141.0) 


91.5(60.0-141.0) 


0.087 


Total testes volume (mL) 


50.0(32.0-70.0) 


50.0(34.0-64.3) 


0.435 


Semen volume (mL) 


3.4(1.3-6.4) 


3.2(1.3-6.5) 


0.781 


Sperm concentration (10^/mL} 


64.3(9.2-225.4) 


62.8(8.5-201 .6) 


0.235 


Total sperm count (1 o'/ejaculate) 


219.5(19.1-780.8) 


197.1(18.6-691.9) 


0.112 



^Data presented as median(5-95^ percentile). 

"^Non-parametric Mann-Whitney two-tailed U-test implemented in PASW software Grad Pack 18.0; asterisk (*) marks P-values resistant to Bonferroni correction for 
multiple testing. 

doi:10.1371/journal.pone.0094244.t003 



with the wild-type homozygote genotype GG-AsnAsn (mean 44.3 
vrs 50.8 mL; P<0.1). 

Among the Estonian oligozoospermic infertility patients, none 
of the four FSHR haplotypes appeared to systematically modulate 
their reproductive physiology (Table 4). We observed single 
genetic associations for the two haplotypes with intermediate 
prevalence (G-Ser, A-Asn), but these cannot be equivocally 
explained by the deviated reproductive physiology in infertile 
men, and the significant P-values (0.01<P<0.05) may have been 
detected by chance only. Infertile patients carrying FSHR G-Ser 
haplotype (frequency 31.2%; high gene expression, but less 
sensitive receptor) had higher serum FSH and lower Inhibin B 
levels, whereas haplotype A-Asn (14.0%; low gene expression, but 
more sensitive receptor) was associated with larger total testes 
volume (Table 4). 

Phenotypic variance of reproductive parameters in 
young nnen explained by FSH-action modulating genetic 
variants 

The contribution of the three FSH-action modulating genetic 
variants {FSHB -21 IG/T, FSHR -29G/A and c.2039A>G) to the 
measured population variance in serum reproductive hormone 
levels, total testicular volume and sperm concentration was 
estimated in the Baltic young male cohort using the REML 
analysis implemented in GCTA software [33]. Here, we took 
advantage of our previously pubhshed datasets on the FSHR 
c.2039A>G (p.Asn680Ser; Table SI; [12]) and the FSHB - 
211G/T (Table S2; [7]) genotype data of the same samples. 
Together, the three SNPs explained 2.3%, 1.4% and 1.0% of the 
phenotypic variance in circulating FSH, serum Inhibin B and total 
testosterone levels, respectively, with the FSHB -21 IG/T account- 
ing for the largest proportion of the variance (serum FSH, 1.4%; 
serum Inhibin B, 0.9%; total testosterone, 0.6%) (Figure ID). For 
the total testes volume, the FSHB and FSHR SNPs together 
accounted for 1.1% of the phenotypic variance. The FSHB and 
FSHR variants almost equally contributed to the variance in total 
testes volume {FSHR c.2039A>G, 0.5%.; FSHB -211G/T, 0.4%; 
FSHR -29G/A, 0.3%). In agreement with the data on individual 
SNPs [7,12], the estimated joint effect on sperm concentration is 
marginal (0.2%), data not shown). 



Discussion 

FSH action in male and female reproductive physiology is 
modulated by genetic variants determining either serum FSH 
levels or the functionality of FSHR [5] . Previously, few small-scale 
studies have inconclusively investigated the effect of the FSHR - 
29G/A alone or in combination with other FSHR genetic variants 
on male quantitative reproductive parameters [18-22]. This the 
first study to show significant associations between the FSHR 29G/ 
A genetic variant and reproductive hormone levels in men. For the 
Baltic young male cohort, we report statistically significant 
association of the FSHR -29 A-allele with higher serum FSH. 
The results fit with the published data on higher gene expression 
level for the FSHR -29 G-allele compared to A-aUele [16], and 
suggests that in the latter case there might be shortage of FSHR 
molecules to bind the circulating FSH. The genetic effect oi FSHR 
-29G/A was enhanced in haplotype-based association analysis, 
which highhghted the FSHR haplotype -29A/2039G (A-Ser; lower 
transcript level and less sensitive receptor) to modulate reproduc- 
tive physiology. It was associated with statistically higher FSH 
(effect 0.40 lU/L), lower Inhibin B (effect -16.57 pg/mL) and 
smaller total testes volume (effect —2.34 mL). Consistendy, the 
most prevalent haplotype, G-Asn, combining the effect of higher 
FSHR gene expression and increased FSH receptor ligand- 
sensitivity was significantly and cumulatively associated with larger 
testes volume (effect 1.34 mL). Thus, our study demonstrates the 
importance of taking into account the cumulative effect of both 
common genetic variants in the FSHR gene, and suggests the 
haplotype-based association analysis instead of single SNP testing. 

To our knowledge, this is the first study aiming at estimating the 
genetic contribution to the normal phenotypic variance in 
circulating FSH level and total testes volume in general 
population. We showed that in healthy young men, the three 
analysed SNPs {FSHR -29G/A, c.2039 A/G; FSHB -211G/T) 
explain together 2.3%, 1.4%, 1.0 and 1.1% of the measured 
variance in serum FSH, Inhibin B, testosterone and total testes 
volume, respectively. Genome-wide association study (GWAS) of 
steroid hormone levels has reported 1 1 SNPs associated with 
estradiol to explain 6.5%, and 6 SNPs associated with testosterone 
to explain 4.4% of the variance in these hormone levels in 
postmenopausal women [34]. A GWAS on male serum testoster- 
one levels reported three SNPs in the SHBG gene to explain 3.8% 
of its variance [35]. 
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A Serum FSH (mean ± SD; lU/L) 



-29 G/A 


G 




A 




l\Sn680Ser 


Asn 


Ser_^^, 


Assn 




Asn 

G 

Ser 


2.9 ± 1.6 
(19.2%) 








3.0+1.6 
(27.6%) 


2.8 ± 1.5 
(9.4%) 






Asn 

A 


3.1 ±1.6 
(14.7%) 


X 


2.9+1.3 

(2.2%) 






3.3 ± 1.8 
(17.3%) 


3.6 ±2.2 
(4.9%) 


3.8 ± 1.8 
(3.0%) 


3.3 ± 1.7 
(1.7%) 


Total testes volume (mean ± SD; mL) 


-29 G/A 


G 




A 




Asn680Ser 










Asn 

G 

Ser 


50.8 ±11,4 
(19.2%) 








48.6 ±10.8 
(27.6%) 


49.2 ± 10.1 
(9.4%) 






Asn 

A 


50.0 ± 9.7 
(14.7%) 


X 


47.6 ±7.9 
(2.2%) 




Ser 


48.9 ± 8.8 
(17.3%) 


47.6 ±10.1 
(4.9%) 


46.3 ±11.8 
(3.0%) 


44.3 ±11.2 
(1.7%) 



Figure 2. Effect of the carrier status of FSHR-29GIIK and Asn680Ser genotype combinations on reproductive parameters. Effect of the 
FSHR -29G/A (rsl 394205) and FSHR Asn680Ser (c.2039A>G, rs6166) genotype combinations on (A) serum FSH level (lU/L; mean ± SD) and (B) total 
testes volume (mL; mean ± SD) among the Baltic male cohort sample (n = 982). The -29G/A and Asn680Ser variants form nine possible FSHR 
genotype combinations: GG-AsnSer (n = 271), GG-AsnAsn (n = 189), AG-AsnSer (n = 170), AG-AsnAsn (n = 144), GG-SerSer (n = 92), AG-SerSer (n = 48), 
AA-AsnSer (n = 29), AA-AsnAsn (n = 22), AA-SerSer (n = 1 7). The prevalence of each genotype combination (%) is shown in brackets. 
doi:1 0.1 371/journal.pone.0094244.g002 



Consistent with the evidence from previous reports in other 
populations [18,19,21,22] and from a recent meta-analysis 
combining seven studies (in total, 1,644 infertility cases and 
1,748 controls) [23], the current study showed no statistical 
difTerence in aUeUc and genotypic distribution oi FSHR -29G/A 
between the Baltic young men cohort (n = 982) and Estonian 
oligozoospermic idiopathic infertility patients (n = 641). Also the 
frequency distribution of the four FSHR haplotypes [formed from 
FSHR -29G/A and c.2039 A/G (p.Asn680Ser)] did not differ 
between the young men and oligozoospermic patients. Thereby, 
we conclude that there is a lack of the association between 
common genetic variants in the FSHR gene and substantial effect 
on male infertility risk. Notably, although the genetic composition 
in young male cohort and oligozoospermic patient group was 
similar, neither FSHR -29G/A alone nor FSHR variants formed 
from the combinations of -29G/A and Asn680Ser appeared to 
signifirandy modulate reproductive physiology among infertility 
patients (Tables 2—4). One possible scenario to explain the 
contrasting outcome of association testing is an important 
hmitation that the hormonal parameters for the two study samples 
had been measured using different methods and laboratories. 
However, hormonal data of both samples were determined in 
accredited centres using standardized commercial assays and have 
been previously successfully applied in genetic association studies, 
providing mutually consistent results [6,8]. A more probable 
scenario supports that in the group of highly selected patients with 
oligozoospermia and male factor infertility, the causes of their 
impaired reproductive function are heterogeneous and alternative 
biological pathways may be involved to maintain their fertility. 
Therefore one should be cautious when interpreting physiologi- 



cally unsound effects of genetic variation on reproductive 
physiology in patients with infertility problems. 

In summary, we conclusively showed the significant effect of the 
FSHR -29 A-allele on male serum FSH level. Effect on FSH and its 
downstream hormonal and testicular parameters were further 
cumulatively modulated by the carrier-status of the FSHR c.2039 
G-allele (p.680Ser). Notably, our study highlights the importance 
in genetic studies of reproductive parameters to utilize population- 
based individuals with normal and undisturbed reproductive 
physiology. We estimated that the three FSH-action modulating 
genetic variants {FSHR -29G/A, c.2039 A/G; FSHB -211G/T) 
account for the substantial proportion of the total normal 
phenotypic variance in male reproductive parameters. Whether 
these polymorphisms may represent genetic risk factors to male 
reproductive disorders apart from fertility, but possibly affected by 
impaired FSH action (e.g. cryptorchidism, hypospadias, testicular 
dysgenesis syndrome and testicular cancer), has to be addressed in 
future studies. 

Supporting Information 

Table SI Marker- trait association analysis and clinical param- 
eters of the two study groups stratified based on the FSHR 
Asn680Ser (rs6166) genotypes of participants. 

(PDF) 

Table S2 Marker-trait association analysis and clinical param- 
eters of the Baltic male cohort sample stratified based on the FSHB 
-211G/T (rsl 0835638) genotypes of the participants. 
(PDF) 
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Table S3 Marker-trait association analysis and clinical param- 
eters of the subgroup of Estonian oligozoospermic idiopathic 
infertility study sample defined according to the World Health 
Organization 2010 criteria (sperm count below 39xl0^/ejaculate; 
n = 408). Parameter data is provided for the subgroups of patients 
stratified based on their FSHR -29G/A, FSHR Asn680Ser and 
FSHB -211G/T genotypes. 
(PDF) 
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